Force localization in contracting cell layers 
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Epithelial cell layers on soft elastic substrates or pillar arrays are commonly used as model systems 
for investigating the role of force in tissue growth, maintenance and repair. Here we show analytically 
that the experimentally observed localization of traction forces to the periphery of the cell layers 
does not necessarily imply increased local cell activity, but follows naturally from the elastic problem 
of a finite-sized contractile layer coupled to an elastic foundation. For homogeneous contractility, 
the force localization is determined by one dimensionless parameter interpolating between linear and 
exponential force profiles for the extreme cases of very soft and very stiff substrates, respectively. 
If contractility is sufficiently increased at the periphery, outward directed displacements can occur 
at intermediate positions, although the edge itself still retracts. We also show that anisotropic 
extracellular stiffness leads to force localization in the stiffer direction, as observed experimentally. 

PACS numbers: 87.10.+e, 87.16.Ln, 87.17.Rt 



Actively generated forces have emerged as a central 
element in the way tissue cells interact with their en- 
vironment, for example during development and tissue 
maintenance [l|. It is becoming increasingly clear that 
cells use physical force to probe the mechanical prop- 
erties of their environments and to create a structurally 
coherent state within the cell and tissue. Experimentally, 
however, it is challenging to measure cellular forces in a 
physiological context. Therefore their effect is usually as- 
sessed in an indirect manner, for example by mechanical 
relaxation after laser cutting Q. 

In order to directly measure the forces that cells ex- 
ert on their environments, two complementary techniques 
have been established over the last decade. Embedding 
marker beads into soft elastic gels enables the measure- 
ment of the displacement field generated by cell traction 
forces applied to the surface of the gel, which can then 
be estimated using elasticity theory [3j. Alternatively, 
microfabrication techniques are used to create an array 
of elastomeric pillars, which for small displacements each 
act as a linear spring 

In order to address the role of forces for larger cell clus- 
ters, both soft elastic substrates and pillar assays have 
been used for confluent layers of epithelial cells. Plac- 
ing such cells on a pillar array, it was found that cellu- 
lar traction forces are localized to the edge of the cell 
layer 0, Q . A qualitatively similar effect has also been 
observed for cell sheets migrating over soft elastic sub- 
strates 0, although in this case the length scales over 
which the stresses decayed were greater. It is tempting 
to assume that force localization to the edge reflects in- 
creased mechanical activity of the cells at the edge, for 
example of putative leader cells 0. However, a more de- 
tailed analysis requires a mechanical analysis of the cell 
sheet coupled to the elastic foundation. 

In this Letter, we analytically solve the problem of a 
contracting cell layer coupled to a layer of springs. For 




FIG. 1. (a) Scanning electron micrographs of an epithelial 
cell layer on an array of microfabricated elastomeric pillars 
(image reproduced from [Eg].) (b) A schematic representation 
of the pillar assay. The cell layer contracts and is resisted by 
a distribution of linear springs on the surface. 

homogeneous contraction, we show that the problem is 
determined by one dimensionless parameter which inter- 
polates between the extreme cases of soft pillars with a 
linear force profile and stiff pillars with an exponential 
force profile. We also show how our results change for in- 
creased contractile activity at the rim and for substrates 
with anisotropic stiffness. 

Elastic model. Microstructured surface assays consist 
of a regular array of flexible elastomeric pillars with de- 
fined surface chemistry, see Fig. [T] For small displace- 
ments and slender pillars, the relation between the trac- 
tion applied to the pillars T and the displacement u fol- 
lows from linear elasticity theory [9( , 

T = feu= (^ p ^)u, (1) 

where r and L are the radius and height of the pillar, 
respectively, and E p is the Young's modulus. The pillars 
thus behave as simple linear springs with effective spring 
constant k. 

In the assays of interest the substrate is usually coated 
with adhesion-promoting molecules so that the cells 
spread thinly and form a confluent layer. We thus as- 
sume the layer behaves as a continuous elastic solid. In 



2 



addition, as the typical lateral extension of the sheet is 
much larger than the layer thickness h, we have the con- 
ditions for plane-stress in the layer and can use a two- 
dimensional model for the cells. The force balance then 
reads 



div • F - fciVu = 0, 



(2) 



where F is the two-dimensional in-plane stress tensor ob- 
tained by averaging the layer stress over the thickness, 
and N is the number density of pillars. For a continuous 
elastic substrate, the restorative force applied to the cell 
layer would not take this simple form. Moreover more de- 
tailed assumptions would be required to define the elastic 
problem of the two coupled elastic layers. Nevertheless, 
the linear traction force term from Eq. [5] can be also 
considered as a first order approximation for soft elastic 
substrates [ToL [Tllj . 

The cell layer is mechanically active, with cells simul- 
taneously repositioning and contracting. Here we con- 
sider only the final configuration of a contractile cell 
layer. Then the problem of modelling layer contraction 
is similar to that of thermoelasticity, where temperature 
changes generate expansion or contraction in the mate- 
rial (lflj. In thermoelasticity one typically introduces an 
expansion term proportional to the target volume change 
into the constitutive relation Q. In analogy we take 
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The first term is the usual plane-stress linear elastic con- 
stitutive relation with E c the cellular Young's modulus, 
v the cellular Poisson's ratio, and the strain tensor. 
The second term in Eq. ([3|) captures cellular contractility, 
where we have assumed that the contraction is isotropic. 
If the layer is under no stress, F^ = and thus P = en 
(summation convention applies), and so P is the local 
target volume change of a material element in the cell 
layer. 

At the layer boundary no stress is applied, giving the 
boundary condition F • n = on T, Integrating Eq. (j5J) 
over the cell area D and applying the zero-stress bound- 
ary condition gives f D kNu — 0. Thus momentum is not 
generated in the system, as expected for a closed system. 

Solution for a contracting stripe. It is instructive to 
first solve our elastic problem in one dimension, although 
experimentally even very anisotropic cell layers always 
contract in both lateral directions. We consider a rect- 
angular sheet with e yy — 0, such that the layer contracts 
only in the x-direction over the stripe width Iq. Substi- 
tuting Eq. (|3J) into Eq. ([2j we find that u(x) satisfies 



d 2 u 1 1. ,dP 



(4) 



with the zero stress boundary condition du/dx = (1 + 
v)Pj1 at x = ±Iq and with u = at x = 0. Here we have 




FIG. 2. (Color online) Plot of (a) displacement \u\ and (b) 
stress F within the cell layer for the stripe geometry with 
Po = 0.7, E c = lOkPa, h = I )im, v = 0.45 and 7 =5, 6, 8, 
10, 15 (from top to bottom in (a) and from bottom to top in 

GO). 



introduced a new length scale, the localization length 
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with the second expression following from Eq. (JT]). 

Assuming that the cell layer contracts uniformly, P 
is constant. For the stripe case, it is convenient to set 
P = — 2Pq/(1+v) where Po is then the target contraction 
in the x-direction. We then find that Eq. §4§ is solved by 



u _ _ p sinh(7x/Z ) 
lo 7COSI17 



(6) 



where we have defined the dimensionless localization pa- 
rameter 7 = lo/l, the ratio of layer dimension to the 
localization length /. The internal stress F(x) can be 
obtained by substituting Eq. ^ into the constitutive re- 
lation Eq. (|3J, see SI 0. 

In Fig. [2] the displacement u and internal stress F are 
plotted for different values of 7. We see that a constant 
cellular contraction leads to non-constant pillar displace- 
ment, with greatest deflections at the edge. The non- 
dimensional parameter 7 controls the profile of the solu- 
tion and the rate of decrease of observed tractions away 
from the layer edge. When 7 is large, the pillars dominate 
the system and restrict the displacements to a small re- 
gion near the rim. When 7 is small, the substrate resists 
less and the tractions decay slower. Note that both dis- 
placement and stress scale linearly with cell contraction 
Po- 

Solution for a contracting disc. We now consider the 
two-dimensional problem of a circular cell layer with 
isotropic contraction, which is directly comparable to ex- 
periments with finite-sized and isotropically contracting 
cell layers. From symmetry arguments u = u(r)e r , where 
r is the radial distance from the centre of the plate. The 
radially symmetric equivalent to Eq. Q reads 



dr 2 



dr 



r \ I, 

F "=2 (1 
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The boundary conditions are now u — at r = 0, while 
at r — r we have P = 2(du/dr + u/r)/(l + v). 

Considering again the case of constant cellular contrac- 
tion, P = —Po, the solution of Eq. ([7} with the specified 
boundary conditions can be calculated in terms of mod- 
ified Bcsscl functions: 



A( 7 ) = (l (7) + ^/i(7) 



(8) 
(9) 



Again 7 = r /l is the ratio of sample dimension and 
localization length, i.e. the decay length scale for the 
traction forces I. From Eq. ([5]) the radial stresses F rr 
and hoop stresses Fee in the layer can be calculated, see 
SI [12 ]. which show a similar behaviour as in the one- 
dimensional case from Fig. EJb), except that the hoop 
stress does not vanish at the boundary. The plot for 
displacement essentially looks like the one for the stripe 
case, Fig. [21[a), see S I 1121 . Our results are consistent with 
the simulations in 13], where a uniformly contracting 
network of actin fibres on a discrete pillar array was also 
observed to exhibit localisation of force to the edges of 
the network. 

Regarding the comparison with experiments, typical 
parameter values are iV = 0.25 /jm" 2 for the pillar den- 
sity, ro = 50 /im for layer dimension, h — 1 /im for layer 
thickness, E c = 10 kPa for cellular Young's modulus and 
v = 0.45 for Poisson's ratio. With pillar stiffness ranging 
over 1 - 200 nN /im" 1 [6(, this gives a typical range for the 
dimensionless parameter 7 ~ 7 — 100. The clear decay 
of the traction profile reported in Q agrees nicely with 
our results. A more spread-out traction profile has been 
reported in Q , which may be attributable to a low value 
of 7. However, as this study focused on the dynamics in 
a migrating sheet, direct comparison is difficult in this 
case. 

To further examine the role of 7 in determining the 
profile of the solution, we consider the limiting cases of 
7 small and large. When 7 is small, the pillars are very 
weak and the layer will be able to contract by a large 
amount. When 7 < 1, we obtain to lowest order a lin- 
ear profile u ~ —P r/2, as one can infer either directly 
from the differential equation or from the Bessel func- 
tion. Alternatively, in the case 7>1, when the springs 
dominate the behaviour of the sheet, the displacement 
localizes strongly to the edge and we get an exponential 
profile, 



P r (l + v) T(1 - r/ro) 
2 7 



(10) 



The exact profile of the solution is thus clearly controlled 
by the non-dimensional parameter, 7, which quantifies 
the relative strengths of pillars versus cells. 

Effect of a contractile rim. Up to now, we have 
assumed that the layer undergoes a uniform contrac- 
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FIG. 3. (Color online) Plot of displacement u in a disc 
with contractile rim for (a) P0/P1 =0.1, 0.3, 0.5, 0.7, 1 with 
n/ro = 0.9 (from top to bottom), (b) n/ro =0.75, 0.8, 0.85, 
0.9, with P0/P1 = 0.5 (from top to bottom). In both plots 
7 = 8, Pi = 0.7, v = 0.45. 



tion. We now apply our model for the contractile cell 
layer to two further situations of large practical inter- 
est, namely non- homogeneous contraction in the layer 
and anisotropic layer stiffness. We first consider what 
the displacements would look like if the layer had a more 
mechanically active rim. We consider that the contrac- 
tion jumps up at the rim: 

p _ f -Pa x < xi or r < n 

\ —Pi X\ < x < Iq or r\ < r < ro ' 

for both the contracting stripe and disc. Eqs. (HI and 
([7]) are unchanged, as are the boundary conditions ap- 
plied at both the origin (u = 0) and edge of the layer 
(zero normal stress). At the interface between the two 
contractile regions (i.e. at x = xi and r — r%, for the 
stripe and disc case, respectively), however, we impose 
the additional conditions that the displacement u and 
stresses F (for stripes) and F rr (for discs) are continu- 
ous across the interface. With these additional boundary 
conditions it is possible to obtain analytical solutions for 
the displacement u, see SI 12|. 

For the disc geometry, the resulting displacement is 
plotted in Fig. [3][a) . Note that the plot of the displace- 
ment u has a kink at the point of transition between the 
two contractilities, and that as the inner region becomes 
weaker and Pq/P\ decreases, this feature is accentuated. 
Equally, as Pq/P\ — > 1, this feature disappears and we 
recover the solution for constant contractility. The dis- 
placements for the contracting stripe solution are quali- 
tatively the same. 

For a sufficiently inactive inner region, there is the 
possibility that contraction of the outer rim dominates 
to such an extent that it can pull the inner region to- 
wards it, resulting in positive displacements at interme- 
diate positions, although the layer as a whole retracts 
its edge. This effect depends sensitively on the spatial 
extent of the contracile rim. For the stripe case with 
xi/xq = 1 — e and e <C 1, positive displacements are 
possible only when P — e 2 7 2 Pi/2, see SI 12]. This de- 
pendence on the extent of the more contractile region 
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FIG. 4. (Color online) Plot of (a) displacement |u|, arrows 
indicate displacement directions and magnitude, (b) traction 
force f/Kro exerted on pillars. In both plots 7 = 5, Po — 0.8 
and v — 0.45, e — 3 (major axis in s-direction), the dotted 
black line indicates the initial shape of the layer. 

is also seen in Fig. G^b). In summary, the existence of 
kinks and positive displacements in experiments would 
be a clear signature of increased contractility at the rim. 

Solution for anisotropic pillars. In a natural setting, 
cells are often confronted with highly anisotropic environ- 
ments, which can have a strong effect on their response. 
This important effect can be investigated experimen tally 
by using pillars with elliptical cross-sections 0, [|1 [l4j . 
Then the apparent stiffness of the substrate depends on 
the direction, 

k{6) = K (1 + ecos 2 6») , (12) 

where 6 is defined to be the angle of the pillar deflection 
with respect to the major axis a. The minor axis is given 
by the relation a 2 = b 2 (l + e). The spring constant is now 
K = 3irE p b 3 a/4:L 3 and so we redefine 7 2 = r%KN(l - 
v 2 )/hE c . Using these arrays it is observed in [y, [l4[ that 
tissues elongate along the major axes of the pillars, with 
cellular force being concentrated at the pointed ends of 
the layers. 

Combining Eqs. © and © with k given by Eq. (1T2"1) . 
and with the boundary conditions as before, the prob- 
lem is now fully two-dimensional. We consider again the 
paradigm problem of an initially circular layer, contract- 
ing isotropically, and solve the system numerically us- 
ing finite element methods for pillars with a/b = 2 (i.e. 
e = 3) and the major axis in ^-direction (MATLAB PDE 
toolbox, The Mathworks, Natick, MA). For isotropic pil- 
lar elasticity, our numerical solutions agree perfectly with 
the analytical results given above. The numerical results 
for the anisotropic case are plotted in Fig. 0] Due to 
the symmetry of the problem, it is sufficient to solve the 
problem on a quadrant, imposing zero x, y displacement 
on y = 0, x = 0, respectively. We see that an initially 
circular layer on the anisotropic surface can be induced 
by uniform isotropic contraction into a shape that is elon- 



gated in the direction of larger stiffness (x-direction) . We 
also see that the traction forces / = k(9)u are concen- 
trated at the pointed ends of the layer. We do not expect 
this effect to be sufficient to account for the shape of the 
elongated tissue islands reported in [(J [l4| , but it is im- 
portant to note that an anisotropy in tissue shape and 
force profile can be induced purely through the substrate 
properties. Indeed this result of uniform cellular tension 
generating anisotropic internal stresses could be used by 
cellular assemblies to sense the direction of maximal stiff- 
ness and to polarize accordingly. 
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